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1.  Introduction 

For  the  most  part,  the  atmosphere  provides  a  reliable  media  for 
signal  propagation  in  radio  communications  systems.  It  is  largely 
homogeneous  and  can  be  modelled  as  free  space  for  most  situations.  When 
signals  pass  through  the  ionosphere,  however,  a  great  deal  of  complica¬ 
tions  arise. 

Radio  waves  undergo  scattering  in  a  structured  plasma,  as  in  the 
ionosphere.  The  resultant  perturbations  imposed  on  the  signal  can  prove 
deleterious  in  a  communications  system.  Temporal  fluctuations  in  the 
received  signal  translate  into  distortion  after  demodulation.  Both 
amplitude  fading  and  phase  fluctuations  affect  various  types  of  systems. 

In  order  to  characterize  these  fluctuations,  continuous  wave  signals 
are  transmitted  by  satellite  *md  received  on  the  ground  after  passing 
through  the  ionosphere.  In  this  report,  data  from  the  DNA  wideband 
satellite  experiment  were  obtained  at  various  stations  around  the  world 
run  by  SRI  International.  The  received  complex  signal  displays  the 
variation  of  the  trans ionospheric  scintillations. 

By  studying  the  temporal  structure  of  the  ampltide  and  phase 
scintillations  insight  can  be  gained  into  the  structure  of  the  irregu¬ 
larities.  The  actual  data  available  are  from  low,  medium  and  high 
latitudes.  At  equatorial  and  auroral  latitudes  there  is  considerable 
ionospheric  activity  contributing  to  extensive  scintillation.  Close 
scrutiny  is  given  to  these  data  sets. 

In  addition  to  interpreting  the  scintillation  data  in  terms  of 
geophysical  parameters  in  the  ionosphere  (such  as  mean  square  deviation 
of  electron  density,  the  irregularity  spectrum,  etc.)  the  statistical 
studies  of  the  scintillation  data  can  also  yield  information  that  is 
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useful  to  communication  engineers.  For  example,  the  high  frequency 
asymptote  of  the  amplitude  power  spectrum  is  found  to  follow  a  power 
law  of  the  form  v  n,  which  predicts  a  frequency  dependence  of  the 
scintillation  index  in  the  form  f  (n+^)/A  in  ty,e  weak  scatter  regime. 

Thus  if  n  is  determined  experimentally,  the  dependence  of  scintillation 
on  radio  frequency  f  is  known  at  least  when  the  scintillation  is  not 
strong.  As  an  another  example,  the  value  of  correlation  time  is  a 
desirable  information  to  have  in  applications  of  time  diversity.  If 
the  temporal  variation  in  the  signal  is  created  by  drifting  of  frozen- 
in  random  field  as  usually  is  the  case,  the  correlation  time  can  be 
easily  translated  to  the  spatial  correlation  interval.  Of  course, 
in  applications  of  spatial  diversity  it  is  desirable  to  know  the  value 
of  spatial  correlation  interval.  The  four  ground  stations  from  which 
data  were  collected  are  Poker  Flat,  Alaska  (auroral  station),  Stanford, 
California  (mid-latitude  station),  and  Ancon,  Peru  and  Kwajalein,  Marshall 
Island  (both  equatorial  stations) .  In  the  natural  environment  there  is 
least  scintillation  at  mid-latitude,  followed  by  the  auroral  station. 

The  equatorial  stations  have  the  strongest  scintillation.  Consequently, 
in  this  sense  the  equatorial  environment  is  the  closest  simulation  to 
the  case  of  artificially  modified  ionosphere. 

Previous  work  in  this  area  has  been  carried  out  with  a  different  type 
of  data.  Umeki  et  al.  (1977a, b)  and  Myers  (1979)  used  data  from  the  ATS-6 
Radio  Beacon  Experiments.  Frequency  coverage  in  this  data  was  limited  to 
AO,  140,  and  360  MHz,  with  the  latter  providing  the  phase  reference.  In 
addxcion,  the  satellite  was  geostationary  and  the  data  analyzed  provided 
temporal  latitude  coverage  only.  Orbiting  satellite  data  at  many  frequencies 
from  diverse  geographic  locations  is  used  here  to  overcome  these  limitations. 
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Previous  studies  of  the  SRI  data  have  also  been  carried  out.  Fremouw 
et  al.  (1977)  undertook  studies  using  large  amounts  of  data  to  draw  general 
conclusions.  In  this  report  specific  data  sets  will  be  examined  extensively 
in  "case  study"  analysis. 

Chapter  2  gives  some  motivation  for  the  subsequent  studies  by  predicting 
the  theoretical  bases  for  scintillation.  Weak  and  strong  scintillation,  in¬ 
cluding  a  general  theory  on  the  autocorrelation  of  the  complex  signal,  are 
examined.  Chapter  3  is  devoted  to  the  techniques  used  in  acquisition  and 
analysis  of  the  data.  Results  and  conclusions  are  then  presented  in 
chapters  4  and  5,  respectively. 


% 
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2.  Scintillation  Theory 

In  this  chapter,  some  of  the  physics  of  a  model  ionosphere  will  be 
discussed  in  order  to  develop  predictions  about  its  behavior. 

2.1  Weak  Scintillation 

For  our  purpose  of  scintillation  study,  the  ionosphere  can  be  thought 
of  as  a  slab  with  thickness  L  containing  electron  density  irregularities 
conducive  to  the  scattering  of  radio  waves.  An  orbiting  satellite's  trans¬ 
missions  are  received  at  an  earth  station,  thus  exhibiting  the  effects  of 
spatial  variations  in  the  slab  (Hg.  2.1.1).  The  electric  field  in  the 
slab  is  given  by  the  solution  of  the  inhomogeneous  wave  equation. 

V2E+k2e(r)E  =  -V(E«V«,n[£  (r)  ] )  (2.1.1) 

k  =  X  =  wavelength 

e  *  dielectric  permittivity 

Because  compared  to  the  velocity  of  the  radio  signal  (as  well  as  the 
satellite)  the  variations  in  the  irregularity  change  very  slowly,  the  per¬ 
mittivity  can  be  considered  a  function  of  position  only.  It  can  be  formulated 
as 

e(r)  =  <£>[!+£!(?)]  (2.1.2) 

where  £j(r)  is  the  random  component  of  the  permittivity. 

For  weak  scintillation,  or  single  scattering,  the  Rytov  approximation 
can  be  applied  (Tatarski,  1971).  For  a  large  distance  between  the  transmitter 
and  the  top  of  the  slab,  plane  wave  incidence  of  the  form  ?oexp(-jic*r)  can  be 
assumed.  This  gives  a  Rytov  solution  to  (2.1.1)  as  (Tatarski,  1971) 


Transmitter 


Figure  2.1.1  Geometry  of  Model. 


E(r)  -  £0exp(-j£*r+\jj(r)) , 


(2.1.3) 


4>(r)  -  rg  HI  ANe(r')exp[j^*(r-r')]exp[jk|rT?' |)(-^5fry) 


r  ■  classical  electron  radius 
e 


AN£  =  electron  density  fluctuation 


The  complex  quantity,  ip(r) ,  can  be  written  as  the  sum  of  its  real  and 


imaginary  components. 


<l>(  r)  =  x(^)+JS(r) 


(2.1.4) 


The  real  part,  x»  will  yield  amplitude  variations  while  the  imaginary  component, 
S,  represents  the  phase  fluctuations. 

It  has  been  shown  (Wernik  et  al.,  1973)  that  at  the  receiver  on  the  ground 


f  4  “ 

<X2>  *  ?  '2  L  |  FF(K)*N(<)icdK, 


FT7(k)  =  1  -  sin  ("2j~)  cos  [—(z  -  T)) 


(2.1.5) 


and  also 


where 


~-f  ~  |  F$(K)$N(K)Kd<, 


V°  "  1  +  sin  (2T}  cos  [r(z  ■  2)3 


ic  •  spatial  wave  number 

$^(ic)  *  power  spectrum  of  electron  density  fluctuations 
f  *  signal  frequency 


(2.1.6) 


f  ■  plasma  frequency  of  background  ionosphere 
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Experimentally  (Rufenach,  1972;  Crane,  1976),  a  power  law  spectrum  has 
been  demonstrated  for  the  irregularities;  i.e.,  4>n<*k~p.  The  specific  form 
used  is 


_ <(AN/Nn)2>nrn3  T(p/2) _ 

2irr(3/2)r(£^-)  [l+r02(K2  +ti2k  2)]p^ 
^  x  y 


(2.1.7) 


<(AN/Nq)2>  *  variance  of  electron  density  fluctuations 
n  =  anisotropy  factor 
rg  =  outer  scale  size  of  irregularities 
F(x)  =  gamma  function 

The  temporal  power  spectra  of  the  received  signal  can  now  be  predicted. 
Assuming  frozen-in  irregularities,  Umeki  et  al.  (1976)  gave  an  amplitude 
spectrum  of 


PA(°) 


irk2L 

4v 


FF(q)*N(K)dicy 


(2.1.8) 


The  phase  spectrum  is  (Singleton,  1974) 


Vn) 


irk2L  ,£2.-  4 


4v 


F$(q)$N(K)dKy 


(2.1.9) 


where  =  2irf 


f 

v 


n2  .  2 

v‘  y 

temporal  frequency 
satellite  velocity, 


projected  onto  slab 


For  amplitude,  the  temporal  spectrum  begins  to  roll  off  at  the  Fresnel 
frequency,  f  “vM...  Given  the  condition  of  power  law  irregularity  spectrum, 

f  r 

the  temporal  amplitude  power  spectrum  will  decrease  as  at  the  high 


frequency  end.  On  the  other  hand,  the  temporal  phase  spectrum  should  be 
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proportional  to  f  ^  ^  for  all  frequencies  above  the  cutoff  used  to  detrend 
the  data.  The  low  frequency  leveling-off  characteristic  of  the  amplitude 
spectrum  is  due  to  the  Fresnel  filtering  effect.  Irregularities  larger  than 
the  Fresnel  size  (X  **/X(2z-L))  do  not  contribute  significantly  to  amplitude 

r 

scintillation. 

2.2  Strong  Scintillation 

The  above  derivations  are  valid  only  under  the  conditions  |ej|«l,  <x2><<1» 
and  <S2>«1.  These  no  longer  hold  for  strong  scintillations.  Multiple  scattering 
effects  present  an  extremely  complicated  situation.  Rumsey  (1975)  and  Marians 
(1975)  have  carried  out  some  theoretical  computations  for  the  amplitude  spectrum 
using  a  phase  screen  model.  Their  results  predict  that,  even  under  conditions 
of  strong  scintillation,  the  power  spectrum  will  roll  off  with  a  constant  slope. 
Liu  et  al.,  (1974)  and  Yeh  et  al.,  (1975)  have  studied  the  saturation  effects 
of  the  amplitude  scintillation,  frequency  dependence  of  scintillation  index 
and  the  correlation  time  under  multiple  scattering  conditions.  Some  of  those 
results  will  be  applied  to  the  interpretation  of  the  data. 

2.3  Correlation  Characteristics 

Some  current  theoretical  investigations  are  predicting  behavior  that 
applies  both  in  strong  and  weak  cases.  In  particular,  both  the  average  field 
and  the  mutual  coherence  function  of  the  field  can  be  obtained  analytically 
under  the  forward  scatter  approximation. 

Let  us  define  the  complex  amplitude  of  the  field,  U,  by  E*Uexp(-jkz) . 

Insight  can  be  gained  into  the  actual  structure  of  the  irregularities  by 
observing  the  autocorrelation  of  the  received  signal,  U(z,p) | z*ground.  The 
mutual  coherence  function  of  the  complex  amplitude  is  defined  as 


r2  *  <U(z,p i)U*(z,p2)>  (2.3.1) 

For  plane  wave  incidence  a  solution  can  be  obtained  for  both  <U>  and 
r2  (Tatarski,  1971)  inside  the  irregularity  slab: 

<U>  =  Uoexp[-  A(0) ]  (2.3.2) 

r2  =  U02exp  {-—  [A(0)-A(p) ] } ,  (2.3.3) 

A(p)  =  2  dimensional  correlation  function 
for  ej  irregularities 

CO 

=  |  Be(p,z)dz 

—CO 

Uq  =  incident  field  amplitude 


Below  the  slab,  z>L,  the  results  are 

r2  ■  Uo2exp{-  -pp  (A(O)-A(p)]}  (2.3.4) 

<U>  =  u0exp[-  -^jpA(O)  ]  (2.3.5) 

Define  L=£n[r2/<U>2]  and  hence 

L  -  X,nr2-Jln<U>2  (2.3.6) 

Substitution  of  (2.3.4)  and  (2.3.5)  in  (2.3.6)  yields 


i  /~y\  k2  L  . 
L(o)  *  — 7—  A(p ) 


(2.3.7) 
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In  terms  of  the  percentage  electron  density  fluctuation  AN/N0 ,  (2.3.7)  can 
be  written  as 

_  ir2f  4L 

L(p)  -  —  •£  2  A^(p),  (2.3.8) 

t 

c  ■  speed  of  light 

Afl  “  two-dimensional  transverse  correlation  function 
of  AN/N0 


Assuming  "frozen  in"  irregularities  (Tatarski,  1971),  the  transverse  dis¬ 
placement  becomes  p  «vt.  Thus,  the  correlation  of  the  irregularities  can  be 

observed  directly  from  the  logarithm  of  the  correlation  of  the  received  complex 
signal. 

In  terms  of  the  irregularity  spectrum,  *N,  A^v.t)  can  be  formulated  as 

A^v.t)  -  2tt  ||  4>N(0,K)e_jK‘vtdK  (2.3.9) 

Next  we  evaluate  the  temporal  spectrum  of  A^v.t) 

j  Vv.t)*^  dt  *  2tt  J|  *n(0,£)  d<  |  ej^-v)t  d|. 

*  ^ 2  ||  <I,N(0,Ic)dK6(fi-K«v)  (2.3.10) 

x>  4tt2l 


Let  “y»  «  "  then  I  becomes 


y  y  y  v  F  y  v 
y  y 


*N(0*VKy)5(!!-KxVy)dKxdKy 


||  *N^0,lcx*  v  ^<5^~icxvx-y)d'cxdy/v 


(2.3.11) 


1  (  n~Kx\ 

-y  j  V°’V  —f-*>  «X 
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Further  evaluation  depends  on  the  behavior  of  $N.  As  an  example,  let  us 
consider  the  case  of  isotropic  irregularities  such  that  (2.1.7)  reduces  to 

(2.3.12) 


*N  *  Cn<1+r02K2) 


-(V+|) 


V  = 


(P-1) 

2 


C..  = 


:(AN/N0)2r03r(P/2) 


"  27rr(3/2)r(-^) 

Furthermore,  let  us  assume,  without  loss  of  generality,  that  all  velocity  is 
in  the  y  direction  (v  =0) .  This  gives 


ft-K  v  rn2fi2 


-(\Hrr) 

$  (0,K  , - 2-2).  r  [1+ - jr—  +  r02<  2]  2 

N’x  v  '  N  v  ^  u  x 

y  y 


(2.3.13) 


and  (2.3.11)  becomes 


T  _  CN  2 

I  -  " —  I  (1+  — -j—  +  rn‘ic  /)  die 

u  x  x 


(2.3.14) 


2  CN  (v-1) ! 


ro 


2n2 


-v 


((2v-n:P (1+  v  2  ) 


'  (2v-l) 

rovy  y 


where  v  is  assumed  to  be  an  integer  and  (2v-l)!!  =  1*3*5  ...  *(2v-l).  From 
(2.3.8),  (2.3.10)  and  (2.3.14)  we  can  compute  the  spectrum  of  L. 

TT2f  kL 

|  dt 

(2.3.15) 


TT2f  4 L  f 

\  =  J  Vv.o 


4TT2f  4  C  v  (v-1)  : 
0  <■ 


fl2rn2  'V 


o'  V(2v-1)!.,)  (1+  ~V^~) 


Thus,  for  isotropic  irregularities,  the  power  spectrum  of  the  quantity 
L  defined  in  (2.3.7)  yields  direct  information  about  the  spectral  shape  of 
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the  density  fluctuations.  For  a  more  general  case,  (2.3.11)  can  be  evaluated. 

These  theoretical  results  for  both  weak  and  strong  scintillations  will 
be  used  later  In  the  thesis  for  purposes  of  data  Interpretation. 


3.  Data  Acquisition  and  Processing 


In  this  section  the  techniques  used  for  acquisition  and  analysis  of  the 
data  will  be  explored. 

3.1  Design  of  Experiment 

In  May  of  1976  Stanford  Research  Institute  initiated  the  DNA  Wideband 
Satellite  Experiment  for  the  Defense  Nuclear  Agency.  Satellite  P76-5  was 
launched  into  a  1000  km  near-polar  orbit.  The  satellite  carried  a  multi¬ 
frequency  coherent  transmitter  radiating  10  CW  signals  ranging  from  VHF  to 
S-band  (Table  3.1.1).  The  program  was  designed  to  study  the  scintillations 
of  radio  waves  induced  by  the  ionosphere  (Fremouw,  et  al.,  1978). 

TABLE  3.1.1  Data  Channel  Frequencies 
Channel  Frequency  MHz  Harmonic  (of  11.4729  MHz) 


VHF 

137.6748 

12 

UL3 

378.6057 

33 

UL2 

390.0786 

34 

UL1 

401.5515 

35 

UC 

413.0244 

36 

UU1 

424.4973 

37 

UU2 

435.9702 

38 

UU3 

447.4431 

39 

L-band 

1239.073 

108 

S-band 

2891.171 

252 

Data  was  collected  at  four  ground  stations.  At  mid-latitudes,  a  receiver 
was  located  at  Stanford,  California.  For  equatorial  observations,  stations 
were  located  at  Ancon,  Peru  (near  Lima)  and  Kwajalein  atoll  in  the  Marshall 
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Islands.  Study  of  auroral  behavior  was  facilitated  by  a  receiver  at  the  Poker 
Flat  Rocket  Range  near  Fairbanks,  Alaska  (Table  3.1.2). 

TABLE  3.1.2  Ground  Station  Locations 


station 

Geodetic 

latitude 

Coordinates 

longitude 

Geomagnetic 

latitude 

Poker  Flat,  Alaska 

65°7'34"N 

147°29*14"W 

65.4°N 

Ancon,  Peru 

11°46’34"S 

77°9 '00"W 

0.4°N 

Kwajalein,  Marshall 
Islands 

9°24'08"N 

167°28'10"E 

4.4°N 

Stanford, 

California 

37°24 ' 10"N 

122°10'27"W 

42.8°N 

The  receivers  are  designed  to  yield  a  signal  at  the  differential  doppler 
frequency  between  the  appropriate  measured  signal  and  the  S-band  reference 
(Fig.  3.1.1).  A  synthesizer,  phase-locked  to  the  reference  phase,  provides 
quadrature  outputs  to  the  measured  signals  phase  detectors.  These  (x,y)  out¬ 
puts  are  bandlimited  to  150  Hz,  digitized,  and  recorded  at  500  samples  per 
second . 

Initial  processing  is  performed  in  the  field.  The  data  is  decimated  to 
100  Hz  samples  and  converted  to  intensity  and  phase.  In  order  to  remove  trends 
due  to  TEC  a  detrending  procedure  is  applied.  First,  the  intensity  and  phase 
are  smoothed  by  a  0.1  Hz  (ten  second)  6-pole  Butterworth  low-pass  filter.  De¬ 
trended  intensity  is  obtained  by  evaluating  the  quotient  of  the  original  signal 
and  its  trend.  The  detrended  phase  results  from  subtracting  the  trend  from 
the  raw  data.  The  data  are  then  scaled  to  give  unity  mean  intensity  and  zero 
mean  phase  (Fig.  3.1.2).  These  signals  are  then  recorded  on  nine-track  tape. 
Associated  summary  tapes  are  also  generated  containing  second  order  signal 
statistics  and  geometry  information  on  the  satellite  pass. 
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After  examining  strip  charts  and  signal  statistics,  several  of  these  tapes 
were  copied  from  the  SRI  data  library.  Each  tape  contains  one  satellite  pass 
of  approximately  fifteen  minutes.  The  ten  second  detrended  tapes  contain  nine 
channels  of  information:  three  files  with  three  channels  each.  The  channel 
frequencies  used,  available  from  each  station,  were  137.7  MHz  (VHF) ,  378.6  MHz 
and  447.5  MHz  (UHF) ,  1239.1  MHz  (L-band)  and  2891.2  MHz  (S-band).  For  each 
pass  copied,  the  corresponding  summary  tape  was  also  copied. 

3.2  Initial  Conversion 

Before  the  data  could  be  useful,  it  had  to  be  converted  to  a  form  access¬ 
ible  to  the  Cyber  175  at  the  University  of  Illinois.  This  conversion  involved 
essentially  two  tasks.  The  numerical  values  had  to  be  unpacked  and  converted 
from  SRI's  (Hewlett-Packard)  formats  to  the  data  representations  appropriate 
for  the  Cyber.  In  addition,  the  file/record  structure  was  organized  into  a 
more  useful  form. 

The  HP  formatted  data  consists  of  two's  compliment  integers  and  double¬ 
precision  floating  point  words  (Fig.  3.2.1).  Each  word  contains  16  data  bits. 
The  conversion  process  maps  these  to  60  bit  one's  compliment  integers  and 
floating  point  Cyber  words  (Fig.  3.2.2). 

The  format  of  the  ten  second  detrended  tapes  obtained  from  SRI  is  cumber¬ 
some.  For  the  computations  in  the  subsequent  analyses  it  is  necessary  to  have 
the  data  files  divided  by  frequency  and  subdivided  into  magnitude  and  phase. 

The  SRI  format  intermixes  the  different  frequencies  (channels)  and  their  re¬ 
spective  magnitudes  and  phases  (Fig.  3.2.3)  within  each  record. 

The  amplitude  (magnitude)  of  the  signal  is  more  useful  for  our  purposes 
than  the  intensity.  This  is  computed  by  taking  the  positive  square  root  of 
the  intensity.  Because  of  the  very  large  number  of  data  points,  they  were 
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Fig.  3.2.1  HP  Numerical  Representations 
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Fig.  3.2.2  Cyber  Numerical  Representations 


records  1  and  2  duplicates  of  raw  data  records 


records  3-end  of  file  (between  1000  and  2000  records): 


word  1-2 
3-4 
5-6 
7-8 
9-10 
11-12 
13-600 
601-720 
721-727 
728-729 
730-742 


'\ 

J 

) 

} 


channel  1  ''N 


channel  2 


channel  3  J 


intensity 
phase 
intensity 
phase 
intensity 
phase 

repeat  (50  times)  above  sequence 
low-pass  reference  signal 
time  information 
antenna  controller  position 
blanks 


^  double-precision  floating 
point  detrend  values 


Fig.  3.2.3.  SRI  Format  for  10-second  detrend  tapes  used 
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split  and  recorded  in  records  of  4000  points,  maximum.  Thus,  an  array  of  m 
points  will  occupy  n  records  where 


m 


4000 


n  = 


^4000 


+  1 


,  mod4000(m)  =  0 


>  mod4000<m)  *  0 


(3.2.1) 


The  reformatted  files,  with  separate  amplitude  and  phase  for  separate 
channels,  are  in  the  form  shown  in  Fig.  3.2.4.  In  any  given  file,  however. 


two  channels  may 

be  omitted 

records  1-n 

amplitude 

n+l-2n 

phase 

2n+l-3n 

amplitude 

3n+l-4n 

phase 

4n+l-5n 

amplitude 

5n+l-6n 

phase 

} 

) 

} 


channel  1  *\ 


channel  2  )  CYBER  real  values 


channel  3  J 


Fig.  3.2.4.  Reformatted  10-second  detrend  file 


The  data  sets  are  named  through  a  simple  procedure.  Each  satellite  is 
referred  to  by  the  first  letter  of  the  receiving  station  and  an  ordinal  index. 
K1  is  the  first  pass  recorded  from  the  Kwajalein  station  and  P3  the  third  from 
Poker  Flat.  If  more  than  one  segment  of  data  is  used  from  a  given  pass,  a 
second  index  indicating  which  record  the  segment  begins  in  is  assigned.  For 
example,  A4R16  starts  in  the  sixteenth  record  of  pass  A4. 

3.3  Analysis  Procedures 

Of  initial  importance  when  examining  the  data  are  several  second  order 
statistics:  scintillation  index  and  phase  standard  deviation.  The  data 
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is  plotted  to  locate  periods  of  scintillation  and  these  statistics  are  used 
to  verify  it  emperically  and  evaluate  its  stationarity .  Their  dependence  on 
frequency  is  also  of  interest. 

The  scintillation  index,  S^,  is  defined  as  the  normalized  standard  devia¬ 
tion  of  the  intensity. 


<I2>-<I>2 

<I>2 


(3.3.1) 


Likewise,  phase  scintillation  is  displayed  by  the  standard  deviation. 


(3.3.2) 


These  are  computed,  for  discrete  data  points,  using  the  standard  arithmetic 
mean  formula. 

n 

<x>  =  —  £  x.  ,  n  values  of  x  (3.3.3) 

k=l 

After  locating  periods  of  scintillation,  segments  of  stationary  data 

must  be  isolated  for  further  processing.  Stationarity  implies,  ideally, 

that  the  above  second  order  statistics  remain  time-invariant.  For  my  purposes 

I  established  the  criteria  that  the  standard  deviation  of  S,  be  less  than  20% 

A 

throughout  the  duration  of  the  segment.  After  selecting  stationary  data  seg¬ 
ments,  two  types  of  correlation  were  computed;  complex  autocorrelation  and 
real  autocovariance.  The  complex  autocorrelation  requires  that  the  data  be 
in  the  form: 


U  *  X+jY,  X  =  A  costp ,  Y  =  A  sin<|>  (3.3.4) 

The  mean,  <U>,  is  then  computed  as  well  as  the  autocorrelation 


Ry(T)  -  <U(t)U*(t+T)> 


(3.3.5) 
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Also  of  interest  is  the  real  autocovariance.  This  is  computed  individually 
for  amplitude  and  phase.  The  mean  is  removed  and  the  autocorrelation  formula 
applied . 


Ra(t)  =  <(A(t)-<A>)(A(t+T)-<A>)> 
=  <A(t)A(t+t)>-<A>2 


(3.3.6) 


The  real  autocovariances  and  the  logarithm  of  the  complex  autocorrelation  are 
then  saved  for  spectral  computations. 

The  computation  of  power  spectra  is  facilitated  by  the  fact  that  for  a 
given  random  variable,  the  autocovariance  and  power  spectrum  form  a  Fourier 
transform  pair.  Thus,  if  a  signal's  auloeovariance  is  known,  the  power  spec¬ 
trum  can  be  realized  by  evaluating  it  in  the  frequency  domain.  The  computational 
situation  is  somewhat  more  complex  so  I  will  outline  the  details  essential 
for  this  realization.  In  an  ideal  case,  the  input  data  and  the  autocorrelation 
would  be  continuous  and  extend  from  negative  infinity  to  positive  infinity  in 
time.  In  the  practical  situation,  the  data  is  sampled  every  ig  seconds.  The 
data  itself  extends  for  t  seconds  while  its  autocovariance  lasts  for  time 


t^.  These  limitations  reduce  the  computations  to  a  statistical  estimation  and 
present  several  sources  of  error. 

Initially,  an  error  is  introduced  by  the  finite  record  length  in  the  time 
domain.  Confidence  limits  for  the  spectral  estimation  must  be  evaluated 


(Blackman  &  Tukey,  1958).  For  record  length  and  autocovariance  length 

t  the  degrees  of  freedom,  k,  can  be  established, 
m 


2t  * 
n 


T  =  T  —t  X 

n  n  3  m 


(3.3.7) 


m 

* 

The  interval  for  90%  confidence  can  then  be  expressed 


interval  - - dB  ,  k>4 

/k^ 


(3.3.8) 


22 


Another  Important  result  of  a  finite  time  domain  record  is  that  the 
spectrum  is  window  distorted.  The  definite  start  and  stop  of  data  is  actually 
multiplication  by  a  rectangular  window  (Oppenheim  and  Schafer,  1975). 

fl  Oft<r 

”  (3.3.9) 

"■  I 

'vP  t<0,t>T 

m 

In  the  frequency  domain,  this  convolves  a  spectrum  of  the  form  sin  x/x  with 
the  desired  spectrum. 

The  large  sidelobes  of  the  sin  x/x  function  contribute  significantly  to 
distortion  of  the  spectrum  being  analyzed.  Hence,  it  is  essential  that  the 
time  domain  data  be  multiplied  by  a  window  with  smaller  sidelobes  prior  to 
transformation.  Ideally,  the  window's  spectrum  would  be  an  impulse  at  DC. 

For  practical  purposes,  the  Hanning  window  provides  a  good  compromise. 

f  0 , 5+0 . 5cos  (— )  0£t<r 

WH(t)  =/  Tm  m  (3.3.10) 

U  t<0,t>rm 

Even  with  small  sidelobes,  any  DC  spectral  component  will  display  them  as 
distortion  in  the  final  computation.  For  this  reason,  it  is  essential  that 
the  mean  be  removed  from  the  time  domain  data  at  some  point  in  the  processing. 

The  range  of  frequencies  to  be  observed  in  the  power  spectra  arises  as 
a  point  of  interest.  The  sampling  theorem  indicates  that  data  in  discrete-time 
are  reliable  only  for  frequencies  below  (2tg)  \  where  Tg  is  the  sampling  period 
(Lathi,  1974).  In  this  case,  Tg»0.01  seconds  so  50  Hz  is  the  highest  useful 
frequency. 

At  the  low  end  of  the  scale,  the  point  at  which  the  data  becomes  useful 
is  determined  by  the  frequency  resolution,  Af^.  Any  variations  in  an  interval 
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from  0  Hz  to  Af^  are  too  slow  to  be  resolved.  For  this  reason,  meaningful 
spectral  information  is  available  only  past  this  interval.  This  resolution 
is  given  by  (Lovelace  et  al.,  1970). 


Af  =  1/t 
r  m 


So  to  yield  two  decades  of  useful  spectrum 


Af  =  0.5  Hz 
r 

x  =2  seconds 
m 


(3.3.11) 


(3.3.12) 


While  greater  resolution  could  be  obtained,  it  would  decrease  the  confidence, 
as  indicated  above. 

In  order  to  keep  the  spectrum  "smooth"  the  discontinuity  at  time  zero 
must  be  removed.  This  sudden  step  would  obviously  introduce  high  frequency 
components.  To  do  this,  the  data  is  replicated  in  reverse  order  at  the  end 
of  the  record. 

The  total  length  of  the  time  series  record  must  be  a  power  of  two  in 
order  to  implement  the  fast  Fourier  transform  algorithm.  In  order  to  "stretch" 
the  data  to  the  appropriate  length,  zeroes  are  appended  to  the  end  of  the 
windowed  autocovariance.  The  resultant  record  input  to  the  transform  is  of 
the  form  shown  in  Fig.  3.3.1. 

The  effects  of  this  also  include  the  addition  of  a  neater  appearance 
to  the  spectral  plot.  Because  the  spectra  can  be  evaluated  only  at  discrete 
frequencies  the  spacing  of  these  frequencies  becomes  important.  Too  great  a 
distance  adds  uncertainty  to  the  space  between  frequencies  or  the  picket  fence 


effect  (Oppenheim  &  Schafer,  1975).  The  frequency  increment  is  given  by 


SfS  •  Hz 


(3.3.13) 


•'■'***  a  Htomm  imm* 


r  ,na*>,m* 


where  N  is  the  number  of  points  in  the  transform.  N-2"1  -2048  points  proves 
to  yield  good  results. 

Thus,  a  typical  spectral  computation  will  procede  as  follows,  t  seconds 

n 

of  data  are  reduced  to  seconds  of  autocovariance  by  subtracting  the  mean 
and  autocorrelating .  The  result  is  then  multiplied  by  the  Hanning  window  and 
input,  along  with  its  mirror  image  and  zeroes,  to  the  FFT.  The  result  is  a 
spectral  estimation  with: 


(1)  90%  confidence  within  -  dB 

/k^I 

(2)  a  useful  range  from  0.5  Hz  to  50  Hz  and  resolution  of  0.5  Hz 


(3)  values  computed  every  ~(0'  oi)~(~2~o48)  =  0•0'!,9  Hz 

Plots  are  generated  of  the  detrended  dat^,  amplitude  and  phase  scintillation 
statistics,  autocorrelation  and  spectra. 

As  related  in  Chapter  2,  a  power  law  behavior  is  expected  in  the  spectra. 
At  the  high  frequency  end  this  translates  to  a  linear  relationship  between 
the  logarithm  of  the  magnitude  (dB  level)  and  the  logarithm  of  the  frequency. 
Thus,  as  part  of  the  plotting  process,  I  have  fit  lines  to  the  curves. 

A  standard  software  package  is  employed  that  utilizes  the  least-squares 
method  for  linear  regression.  The  plot  is  first  inspected  to  determine  the 
upper  and  lower  frequencies  between  which  the  line  is  drawn.  The  lower  fre¬ 
quency  is  at  the  corner  defined  by  either  the  spectral  resolution,  detrend 
frequency,  or  Fresnel  frequency  in  the  case  of  weak  amplitude  scintillation. 
The  upper  frequency  is  that  at  which  noise  becomes  dominant. 
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4.  Experimental  Results 
4.1  Pass  Characterization 

After  examination  of  the  time  variation  of  second-order  statistics  of  the 
data,  several  relatively  stationary  data  segments  were  selected  from  each 
receiving  station.  The  plots  of  one  such  segment,  ten  second  detrended,  are 
available  in  Appendix  1.  Preliminary  analysis  shows  some  general  attributes 
of  different  geographic  locations  of  reception. 

Scintillation  strength  varies  greatly  at  different  latitudes.  The  equa¬ 
torial  stations  recorded  the  most  acitivity  with  the  Ancon  data  being  the 
stronger  of  the  two.  The  auroral  station  showed  much  weaker  amplitude  scin¬ 
tillation  but  had  phase  fluctuations  on  the  same  order.  The  data  taken  at 
mid-latitude,  however,  showed  very  little  activity. 

At  any  latitude,  scintillation  level  varies  with  frequency.  The  scin¬ 
tillation  index  is  employed  as  an  indication  of  the  level  of  amplitude  fluc¬ 
tuation.  Recall  3.3.1 


?  _  <i2>-<i>; 
1  =  <i>2 


(4.1.1) 


Averaged  values  of  are  plotted  for  each  receiving  station  in  Fig.  4.1.1. 

Some  type  of  f  n  variation  is  observed  immediately.  The  Poker  Flat 
case  shows  an  almost  constant-slope  with  n;  1.29.  Examination  of  the  Ancon 
and  Kwajalein  cases  reveals  that  the  slope  is  actually  changing.  In  the 
limiting  case  it  appears  that  at  saturation  (S^->1)  the  slope  approaches  zero. 
With  very  small  S^,  however,  n-^vl.5.  The  erradic  behavior  of  the  Stanford 
case  is  most  likely  attributable  to  a  noise  level  on  the  same  order  as  that 
of  the  scintillation. 

In  a  similar  manner,  the  frequency  dependence  of  phase  scintillations 


can  be  analyzed.  The  standard  deviation  (3.3.2)  is  used 
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o  *  (4.1.2) 

$ 

Again,  averaged  values  for  each  station  are  plotted  (Fig.  4.1.2). 

In  most  cases,  o,  varies  as  1/f.  Again,  however,  some  levelling 

4> 

off  is  apparent  for  very  intense  activity.  For  the  Stanford  data,  phase 
scintillation  is  much  higher  than  the  noise  level.  This  level  is  easily 
observed  because  the  S-based  phase,  supposedly  identically  zero,  has  a 
standard  deviation  of  a  few  thousandths. 


4.2  The  Average  Field 

The  levels  of  electron  density  fluctuation  can  also  be  observed  for 
various  passes  by  noting  the  mean  of  the  complex  signal  as  in  (3.3.4). 

Recall  (2.3.5): 

i<2t 

<U>  *  UoGxpt - g^A(O)]  (4.2.1) 


This  leads  to 


£n<U>2 


£n{U02exp[-  A(0)]} 

k2L 

2£nU0-  A(0) 


(4.2.2) 


By  observing  from  (2.3.1)  and  (2.3.3) 

r2  =  <U(z,Pi)U*(z,p?)>  -  Uo2exp{-  ~j~~  [A(0)-A(p)]J  (4.2.3) 

and  evaluating  at  p*0 


r2 


I  *  <u*u>  =  u02 

10 


(4.2.4) 


Recalling  that  <U*U>sl  rewrite  (4.2.2)  as 


PHASE  SCINTILLATION 


Phase  Scintillation 
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*n<U>2  -  -  A(0)  (4.2.5) 

from  the  results  of  (2.3.7)  and  (2.3.8)  this  can  be  rewritten  as 

TT2f 

«,n<U>2  = - c^2  ^(0)  (4.2.6) 

Thus,  the  correlation  function  of  the  electron  density  fluctuation,  A^,  is 
available. 

For  a  given  value  of  A^(0)  equation  (4.2.6)  implies  that 

2.n<U>2  “  -L- 

or 

£n<U1>2  f 2  2 
i.n<U2>2  =  (f7) 

This  hypothesis  was  tested  with  several  sets  of  data.  Some  of  the  characteris¬ 
tics  of  these  data  sets  are  given  in  Table  4.2.1. 


(4.2.7) 


(4.2.8) 


TABLE  4.2.1  Mean  Fields  and  Scintillation  Indices 


Pass 

Frequency  (MHz) 

h 

<U> 

£n<U>2 

K2R4 

378 

0.655 

.598-. 081 j 

-1.0283 

447 

0.605 

•772-.025J 

-0.6794 

1239 

0.163 

.969-. 011 j 

-0.063 

K4 

378 

0.549 

.699-. 044 j 

-0.7162 

447 

0.416 

.770-. 042 j 

-0.5227 

1239 

0.121 

.974-. 007 j 

-0.05269 

P3 

378 

0.068 

.862-. 013 j 

-0.297 

447 

0.057 

.899-. 008 j 

-0.213 

1239 

0.018 

.990-. 002 j 

-0.0201 

P3R16 

378 

0.289 

.503-. 08j 

-1.3743 

447 

0.228 

.621-.065j 

-0.9528 

1239 

0.055 

.96-. 017 j 

-0.0816 
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For  each  case  equation  4.2.8  is  applied  yielding  the  results  in  Table  4.2.2. 


TABLE  4.2.2  Ratios  from  equation  4.2.8 


Pass 

(in<U378>2, 

*n<U447>2 

^u447>2 
( - - — — ) 

£n^U1239>2 

K2R4 

1.51 

10.78 

K4 

1.37 

9.92 

P3 

1.39 

10.60 

P3R16 

1.44 

11.67 

expected 

1.40 

7.68 

For  the  low  frequency  case  an  average  error  of  about  3.5%  is  observed.  In 
the  case  of  the  higher  frequency,  however,  about  40%  error  occurs.  This  error 
is  probably  attributable  to  the  fact  that  <U>  is  so  close  to  unity  in  the  L-band 
its  logarithm  is  very  close  to  zero.  This  makes  the  calculations  highly  sensi¬ 
tive  to  even  a  small  error. 

The  next  step  is  to  isolate  the  frequency  independent  information  from  these 
data.  Rewrite  4.2.6  as 

AN(0)fpo4L  =  -  £n<U>2  (4.2.9) 

Now,  for  a  given  pass,  AN(0)fp^4  L  should  be  constant.  Examine  table  4.1.3. 

Especially  for  the  upper  two  frequencies,  the  computed  values  remain  very 
close.  In  addition  comparing  the  computed  value  for  each  pass  with  the  level 
of  scintillation  indicated  in  Table  4.1.1  shows  another  trend.  As  is  to  be 
expected.  The  amount  of  scintillation  varies  with  the  amount  of  electron 
density  fluctuation. 
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TABLE  A. 2. 3.  Density  fluctuation  coefficient 


pass 

frequency  (MH2) 

V°>fP.‘L 

K2R4 

378 

33 

1.34x10 

447 

1.238xl033 

1239 

8.819xl032 

mean 

1.153xl033 

K4 

378 

32 

9.33x10 

447 

32 

9.52xlOJ 

1239 

32 

7.38x10 

mean 

8.74xl032 

P3 

378 

32 

3.87x10 

447 

32 

3.88x10 

1239 

32 

2.8x10 

mean 

32 

3.52xl0J 

P3R16 

378 

33 

1. 79x10’’ J 

447 

33 

1.74x10 

1239 

33 

1.14xlOJ 

mean 

1. 56xl033 
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4.3  Correlation  Characteristics 

The  autocorrelation  of  the  complex  signal  is  displayed  in  Fig.  4.3.1 
for  pass  K4  at  VHF.  It  is  important  to  note  that,  ideally,  the  autocorrela¬ 
tion  is  purely  real.  Comparison  of  the  magnitudes  of  the  real  and  imaginary 
parts  helps  to  indicate  how  "good"  the  data  is. 

Because,  as  related  earlier,  <U*U>  is  set  at  unity  in  the  field  the 
following  result  is  implied. 

{.nr2  |  «  2n<U*U>  5  0  (4.3.1) 

This  prevents  calculations  like  the  ones  carried  out  for  the  average  field. 

Seperate  correlations  for  amplitude  and  phase  raise  some  interesting 
points,  though.  Define  the  correlation  interval  t  to  be  the  time  it  takes 
for  the  autocovariance  to  drop  to  one  half.  Thus,  smaller  t  corresponds  to 
quicker  decorrelation. 

For  amplitude,  a  dependence  of  on  frequency  and  is  noted.  The 
various  passes  were  separated  into  three  groups  based  on  in  the  L-band. 

At  each  frequency,  within  each  group,  the  mean  ta  is  evaluated.  The  re¬ 
sults  are  plotted  in  Fig.  4.3.2.  For  each  group  the  values  are  normalized 
with  respect  to  the  L-band  value. 

The  best  behaved  curve,  that  of  group  3,  is  very  nearly  proportional 
-1/2 

to  f  .  This  is  attributable  to  the  Fresnel  filtering  effect  related  in 
Chapter  2.  For  weak  scintillation  this  inverse  square  relation  is  to  be 
expected  (Umeki,  et  al. ,  1977a). 

For  the  more  intense  cases,  however,  a  marked  decorrelation  at  low 
frequencies  occurs.  The  exact  behavior  is  difficult  to  assess  because  of 
the  gap  between  447  MHz  and  1239  MHz.  In  both  cases  it  seems  that 


NORMALIZED  AMPLITUDE 
CORRELATION  INTERVAL 


Fig.  4.3.2  Amplitude  Decorrelation 
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increases  with  frequency  up  to  some  point  in  the  UHF  band  and  then  decays 
slightly  into  the  L-band. 

A  useful  task  might  be  to  determine  the  value  under  strong  scintilla¬ 
tion  conditions  for  which  a  signal  is  best  correlated.  A  reasonable  hypothesis 
seems  to  be  that  maximum  correlation  occurs  for  signals  with  S^~0.7.  For  group 
two  this  occurs  around  447  MHz  and  for  group  one  between  447  MHz  and  1.2  GHz. 

The  frequency  dependence  of  phase  autocovariance  is  an  entirely  different 
situation.  The  mean  correlation  interval  <t  >,  is  plotted  in  Fig.  4.3.3  for  each 
of  the  three  groups.  The  Fresnel  filtering  is  absent  from  phase  scintillation. 
Thus,  in  the  weak  case  of  group  3,  the  correlation  interval  is  almost  frequency 
independent.  Multiple  scattering  is  evident  in  the  stronger  groups,  however, 
inducing  increased  decorrelation  at  the  low  frequency  end  of  the  scale. 

4.4  Power  Spectra 

Three  types  of  spectra  are  available  at  each  frequency  for  each  data  set. 
These  are  the  power  spectra  of  amplitude,  phase  and  L=Jlnr2.  These  are  notated 
as  P^,  P^,  and  P^.  The  various  spectra  will  be  presented  in  three  ways.  First, 
those  for  an  individual  data  set  from  each  receiving  station  will  be  discussed. 
Next,  different  segments  from  the  same  pass  will  be  examined.  Finally,  several 
spectra  will  be  compared  in  order  to  note  any  geographic  trends. 

Amplitude  spectra  from  pass  K4  are  presented  in  Fig.  4.4.1.  The  primary 
difference  between  the  VHF  and  UHF  channels  involves  the  corner  frequency 
around  1  Hz.  Notice  the  broadening  in  the  VHF  spectrum  as  opposed  to  the  well 
defined  corner  at  UHF.  It  can  also  be  noted  that  the  broader  spectrum  is  asso¬ 
ciated  with  the  higher  scintillation  index. 

According  to  the  results  in  Chapter  2  the  phase  spectrum  should  display 
the  same  asymptotic  slope  as  its  amplitude  counterpart.  Fig.  4.4.2  does  not 
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Fig.  4.4.1  Pass  K4  Amplitude  Spectra 


bear  this  out  when  compared  to  the  previous  plots.  With  scintillation  level 
approaching  saturation,  the  VHF  phase  spectrum  becomes  flatter.  At  UHF,  how¬ 
ever,  the  amplitude  and  phase  slopes  are  much  closer. 

Figures  4.4.3  and  4.4.4  show  the  spectra  at  VHF  and  UHF.  As  predicted 

in  theory,  the  two  are  very  similar  with  the  exception  of  more  noise  at  VHF. 
This  extra  noise  can  be  attributed  to  the  fact  that  the  imaginary  part  of  the 
autocorrelation  is  slightly  greater  at  VHF.  Nonetheless,  comparison  with  the 
UHF  phase  spectrum  demonstrates  a  great  deal  of  similarity.  Therefore,  where 
the  phase  spectrum  is  not  flattened  by  very  intense  scintillation  it  coincides 
well  with  the  power  spectrum  of  the  electron  density  fluctuations. 

The  amplitude  spectra  for  pass  A4R13  are  shown  in  Fig.  4.4.5.  Again,  the 
VHF  spectrum  appears  broader.  Note,  though,  that  the  curve  actually  turns  up 
to  define  a  corner  frequency  whereas  the  UHF  case,  corresponding  to  a  higher 
S^,  has  a  more  obscure  corner.  In  the  L-band,  the  sharpest  corner  corresponds 
to  the  weakest  scintillation.  Differences  in  the  leveling  off  at  the  right 
end  of  the  spectra  are  due  to  differences  in  the  relative  level  of  quantiza¬ 
tion  noise  in  the  data. 

A  dominant  characteristic  of  the  phase  spectra  (Fig.  4.4.6)  is  some  type 
of  ringing.  This  is  most  evident  in  the  L-band.  As  in  the  Kwajalein  case 
strong  scintillation  corresponds  to  a  flatter  spectrum  for  phase  than  for 
amplitude.  In  this  case  it  is  true  for  both  VHF  and  UHF. 

The  spectra  from  pass  P3R16  confirm  the  same  trends.  The  amplitude 
spectra  (Fig.  4.4.7)  show  the  same  broadening  for  more  intense  scintillation. 
While  noise  causes  the  L-band  slope  to  decrease  rapidly,  the  corner  frequency 
becomes  more  clearly  defined  as  scintillation  becomes  weaker.  Also  of  interest 
is  the  fact  that  the  VHF  and  UHF  asymptotes  have  approximately  the  same  slope. 
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The  phase  spectra  (Fig.  4.4.8)  behave  per  expectation.  The  VHF  spectrum, 
a  strong  multiple  scatter  case,  is  somewhat  flatter  than  the  others.  At  UHF 
the  slope  is  in  the  same  range  as  that  of  the  UHF  amplitude  spectrum.  A  confi¬ 
dence  interval  of  about  3  dB  can  help  explain  some  of  the  error. 

VHF  spectra  for  pass  S1R5  show  some  difficulties  typical  of  the  mid-latitude 
data  (Fig.  4.3.9).  Noise,  uncertainty  and  even  some  interference  around  15  Hz 
are  prominent.  Notice,  though,  the  very  sharp  corner  in  the  amplitude  spectrum, 
based  on  earlier  observations,  this  is  appropriate  for  such  a  weak  case.  As 
predicted  for  weak  scintillation,  the  amplitude  and  phase  slopes  are  about  the 
same. 

One  of  the  advantages  of  low  level  activity  is  that  the  data  yields  a 
very  nearly  real  autocorrelation.  The  spectrum  for  VHF  is  displayed  in 
Fig.  4.4.10.  Note  the  strong  similarity  to  the  phase  spectrum.  The  curve 
drops  off  right  at  the  half  Herz  resolution  point  with  about  the  same  slope. 

Some  of  the  small  scale  detail,  including  the  interference,  is  about  the  same. 

Comparison  of  computations  for  different  segments  from  the  same  satellite 
pass  strengthens  some  of  the  previous  conclusions.  Spectra  from  pass  P3  are 
shown  in  Fig.  4.4.11.  Compare  the  VHF  amplitude  spectrum  to  that  of  P3R16 
(Fig.  4.4.7).  The  slopes  are  about  the  same  but  in  the  weak  case,  P3,  the 
corner  frequency  is  more  clearly  defined.  On  the  other  hand,  the  phase  spectra 
are  about  the  same. 

VHF  spectra  for  pass  A4R16  (Fig.  4.4.12)  show  similar  behavior  for  an 
equatorial  pass.  From  the  same  pass,  the  A4R13  VHF  amplitude  spectrum  (Fig.  4.4.5) 
is  not  as  broad.  Correspondingly,  the  scintillation  is  not  as  strong  in  this 
case.  The  computed  slopes  are  not  the  same  but  this  is  because  the  curves  are 
not  linear  for  the  same  range  of  frequencies.  Again,  for  phase  (Fig.  4.4.6) 
the  spectra  are  about  the  same. 
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Fig.  4.4.8  Pass  P3R16  Phase  Spectra 
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The  spectra  for  pass  A4R16  can  provide  another  interesting  comparison 
with  those  of  pass  K2R4  (Fig.  4.4.13).  Both  scintillation  index  and  standard 
deviation  of  phase  are  the  same  for  these  two  passes  which  can  serve  to  illumi¬ 
nate  any  possible  geographic  differences.  The  phase  spectra  are  about  the 
same.  Those  of  the  amplitude  are  slightly  different.  That  for  pass  A4R16  is 
a  little  broader  and  steeper  in  the  high  frequency  end.  Amplitude  spectra  for 
another  pair  of  passes,  K4  (Fig.  4.4.1)  and  A4R13  (Fig.  4.4.5),  are  also  very 
close  with  equal  slopes.  In  this  case  the  Ancon  spectrum  has  the  more  clearly 
defined  corner  frequency.  Because  the  similarity  extends  into  the  higher  fre¬ 
quencies  as  well,  there  seems  to  be  little  spectral  difference  in  receptions 
at  the  two  equatorial  stations  for  data  sets  of  the  same  scintillation  levels. 

While  no  mid-latitude  activity  is  strong  enough  to  compare  to  that  of 
another  station,  some  differences  between  the  auroral  and  equatorial  spectra 
require  examination.  With  the  same  level  of  scintillation  at  VHF ,  the  amplitude 
power  spectrum  for  pass  A3R5  (Fig.  4.4.14)  is  much  broader  than  that  of  pass 
P3R16  (Fig.  4.4.7).  The  phase  spectra,  on  the  other  hand,  are  about  the  same. 

One  might  acknowledge  this  as  a  geographic  difference  but  the  fact  is 
that  no  comparison  can  really  be  made.  Recall  the  results  in  Fig.  4.1.1. 

Because  the  Ancon  scintillation  is  well  into  saturation  at  VHF  the  scintilla¬ 
tion  index  is  not  an  accurate  indicator  of  the  ionospheric  turbulence.  In 
fact,  while  the  auroral  indices  decrease  with  increasing  frequency,  the  S^, 
values  for  the  Ancon  station  remain  constant  or  even  higher  at  UHF  than  VHF. 
Thus,  the  scintillation  environment  at  Ancon  was  far  more  severe  than  that 
at  Poker  Flat  for  these  two  occasions  and  no  conclusion  can  be  drawn  with 
respect  to  geographical  spectral  trends  from  these  two  sets  of  data.  No 
other  sets  with  similar  scintillation  levels  from  the  two  locations  were 
found  in  the  data  for  further  studies  of  geographical  trends. 
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Fig.  4.4.14  Pass  A3R5  Power  Spectra 
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5.  Summary  and  Conclusions 

The  data  from  the  DNA  Wideband  Satellite  Experiment  provide  a  good  basis 
for  the  study  of  transionospheric  scintillation.  Though  stationary  segements 
of  data  are  difficult  to  locate,  the  results  they  provide  are  fairly  consis¬ 
tent.  The  frequency  dependence  of  scintillation,  relationship  of  average  field 
to  electron  density  fluctuations,  characteristics  of  the  autocorrelation  and 
spectral  as  well  as  geographic  trends  are  all  made  available. 

Scintillation  was  shown  to  generally  increase  at  lower  frequency.  Scintilla¬ 
tion  index  varied  as  f~n  with  n  approaching  zero  for  satura  '^n  (S^=l) .  The 
standard  deviation  of  the  phase  varied  as  f  1  but  levels  off  for  the  strongest 
scintillation  cases. 

The  logarithm  of  the  squared  mean  field,  £n<U>2,  was  demonstrated  to  be 
_2 

proportional  to  f  .  This  trend  was  predicted  theoretically  in  equation  A. 1.8. 
Thus,  as  expected  for  the  same  patch  of  ionosphere,  the  data  at  different  fre¬ 
quencies  yield  a  constant  A^O)  (defined  in  equation  2.3.8).  Comparing  different 
values  of  A^O)  with  corresponding  values  verifies  that  stronger  scintillation 
occurs  when  there  is  greater  electron  density  fluctuation. 

The  correlation  interval,  defined  as  the  amount  of  time  the  autocovariance 
remains  greater  than  half,  depends  both  on  frequency  and  scintillation  leve] . 

By  fixing  the  different  dependences  on  frequency  were  observed.  These  can 
be  explained  by  the  Fresnel  filtering  and  decorrelation  due  to  multiple  scat¬ 
tering.  The  phase  correlation  interval  remained  nearly  constant  for  weak  scin¬ 
tillation  but  decreased  in  the  presence  of  multiple  scattering. 

The  effects  of  strong  scintillation  are  evident  in  the  spectral  plots  as 
well.  While  the  amplitude  spectra  are  broadened  for  very  intense  activity, 
the  phase  spectra  become  somewhat  flatter.  In  the  weak  cases,  as  predicted 
theoretically,  the  amplitude  and  phase  have  the  same  spectral  slope. 


As  predicted  theoretically,  the  spectrum  should  indicate  the  spectrum 
of  electron  density  fluctuations.  Its  independence  of  carrier  frequency  as 
well  as  its  similarity  to  the  phase  spectrum  have  been  demonstrated.  Since 
the  phase  spectrum  has  been  used  to  yield  information  about  the  irregularity 
(Rino,  1980),  the  technique  proves  consistent. 

Certain  geographic  trends  were  also  shown.  Scintillation  is  weak  at 
temperate  latitude,  stronger  at  high  latitude  and  extremely  strong  in  the 
equatorial  region.  Scintillation  index  is  highest  at  Ancon  and  Kwajalein 
though  the  standard  deviation  of  the  phase  is  comparable  at  the  Poker  Flat 
station.  For  both  amplitude  and  pha^e  scintillation,  however,  the  Stanford 
data  shows  the  weakest  activity  by  far. 

In  summary,  the  data  confirms  many  of  the  theoretical  predictions.  While 
orbiting  satellite  data  does  not  always  exhibit  stationarity ,  what  is  used  here 
is  very  valuable  due  to  the  large  number  of  frequencies  available,  the  high 
sampling  rate,  and  the  S-band  phase  reference.  Analysis  of  the  results  in 
light  of  the  various  theoretical  models  provides  insight  both  into  the  char¬ 
acter  of  the  scintillation  as  well  as  the  ionospheric  irregularities  themselves. 
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Appendix 

In  this  appendix  data  from  pass  P3R16  are  presented.  The  amplitude  and 
phase  are  both  ten  second  detrended  with  mean  set  at  one  for  amplitude  and  zero 
for  phase.  In  addition,  continuity  has  been  restored  to  the  phase  in  order  to 
avoid  any  2tt  ambiguity.  Some  information  about  the  pass  is  given  in  Table  A1 
and  plots  at  each  frequency  are  presented  in  Figs.  A1-A5. 

TABLE  A1  Pass  Information 


Pass: 

P3R16 

Location: 

Poker 

Flat,  Alaska 

Time: 

day  67 

,  1978  18:37 

:10.06  UT 

Duration: 

40  seconds 

Frequency 

(MHz) 

S4 

a, (rad. ) 
<P 

137 

0.917 

3.33 

378 

0.289 

1.11 

447 

0.228 

0.940 

1239 

0.055 

0.285 

2891 


0.029 


0.002 
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Fig.  A3  0.1  Hz  Detrended  UHF  signal 
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